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ABSTRACT 

Context. A strong toroidal field can exist in form of a magnetic layer in the overshoot region below the solar convection 
zone. This motivates a more detailed study of the magnetic buoyancy instability with rotation. 

Aims. We calculate the a effect due to helical motions caused by a disintegrating magnetic layer in a rotating density- 
stratified system with angular velocity Q making an angle 6 with the vertical. We also study the dependence of the a 
effect on 6 and the strength of the initial magnetic field. 

Methods. We carry out three-dimensional hydromagnetic simulations in Cartesian geometry. A turbulent EMF due to 
the correlations of the small scale velocity and magnetic field is generated. We use the test-field method to calculate 
the transport coefficients of the inhomogeneous turbulence produced by the layer. 

Results. We show that the growth rate of the instability and the twist of the magnetic field vary monotonically with 
the ratio of thermal conductivity to magnetic diffusivity. The resulting a effect is inhomogeneous and increases with 
the strength of the initial magnetic field. It is thus an example of an "anti-quenched" a effect. The a effect is nonlocal, 
requiring around 8-16 Fourier modes to reconstruct the actual EMF based on the actual mean field. 
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1- Introduction 



The magnetic helds in many astrophysical bodies have their 
, origin in some kind of turbulent dynamo. This means that 
• a part of the kinetic energy of the turbulent motions is di- 
' verted to enhancing and maintaining a magnetic field. This 
magnetic field is generally also random, but under certain 
conditions a large-scale magnetic field can also emerge. Here 
by large-scale we mean length scales larger than the energy 
containing scale of the fluid. This can be the case when 
the turbulence is helical, e.g., owing to the simultaneous 
presence of rotation and stratification. 

The evolution of the large-scale magnetic field can be 
described using averaged evolution equations. In the pro- 
cess of averaging, new terms emerge (e.g., the a effect 
and turbulent diffusion) that result from correlations be- 
tween small-scale velocity and magnetic fields. Here one 
usually considers the case where the magnetic fluctua- 
' tions are caused by the fluctuating velocity acting on the 
mean field. However, under certain conditions it might well 
be the other way around. Imagine, for example, the case 
where initially no velocity is present, but there is instead 
a strong large-scale magnetic field the presence of which 
makes the initial state of zero velocity unstable. In that case 
the magnetic field would be responsible for driving veloc- 
ity and magnetic fluctuations at the same time. This type 
of scenario was first simulated in the context of accretion 
discs where the magneto-rotatio nal instability drives turbu- 
lence (|Brandenburg et al] . ll995l ). and later in t he context of 
the m agnetic buoyancy instability with shear (jCline et al.l . 
[2003), which might apply to the overshoo t layer of the Sun. 
It had already been proposed by iMoffattj ()1978D that, once 
the dynamo-generated magnetic field in this layer reaches 



appreciable strengths, the magnetic buoyancy instability 
can set in and govern the dynamics thereafter. The linear 
phase of this instability in a localized flux layer wit h strati- 
fication and rotation was later studied in detail bv lSchmit^ 
(I1984L Il985h . A necessary but not sufficient condition for 
this instability is 



d_ 

dz 



log 



<0, 



(1) 



which essentially means that the magnetic field modu- 
lus B decreases faster w ith height z than the density p. 
iBrandenburg fc Schmitd (|T998 ) performed numerical calcu- 
lations in presence of rotation and determined the a effect 
of the resulting turbulence by imposing an external mag- 
netic field. 

This type of magnetic buoyancy instability is also re- 
lated to the undul atory instab ility in the absence of both 
rotation a ,nd shear (iFan . 20011) and the double diffusive in- 
stability (jSilvers et al.l 1200 9') in presence of shear and no 
rotation. While the focus of the first study has been on the 
formation of flux tubes from a pre-existing toroidal mag- 
netic layer in a stably stratified atmosphere, in the latter 
a magnetic layer was generated from an initially vertical 
magnetic field in presence of strong shear. It was further 
shown that, when the ratio of magnetic to thermal diffusiv- 
ities is sufficiently low, magnetic buoyancy can still operate 
in the tachocline. 

The focus of this work is twofold. Firstly, we want to 
study the nature of the instability at short times, i.e., in its 
initial linear stage. In particular, its dependence on various 
parameters such as magnetic and thermal Prandtl num- 
bers, angular velocity, strength of the initial field, etc, and 
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compare against the linear theory and previous numerical 
work. It can be argued that in presence of rotation this 
instability produces magnetostrophic waves due to balance 
between Coriolis and Lorentz forces. An important result 
highlighted later is that rotation is not vital to the growth 
of this instability. Secondly, we want to study whether this 
instability constitutes a viable dynamo process, so we want 
to measure the mean-field transport coefficients, namely 
the tensors a an d t? using the quasi-kinemai^ ic test-field 
(QKTF) method (ISchrinner etldl. I2005L 12007^ How ever, 
with one exception ( Vermersch fc Brandenbura . |2009() . the 
QKTF has never been applied to the calculation of trans- 
port coefficients in an inhomogeneous turbulence induced 
by the mean magnetic field itself. Therefore we aim to first 
verify the applicability of the QKTF method to this prob- 
lem. For a review on transport coefficients and thei r deter - 
mination using test fields; see iBrandenburg et al.l ()2010D . 
The applicability of this method to problems with an ini- 
tial magn etic field and fluctuations generated from it is dis- 
cussed in iRheinhardt &: Brandenbura (|2010D . 



2. The Model 

We consider a setup simila r to that described in 
IBrandenburg &: SchmittI (|l998f l. The computational do- 
main is a cuboid with constant gravity, gz, pointing in the 
negative z direction, and rotation ft making an angle 9 
with the vertical. The box may be thought to be placed 
at a colatitude 6 on the surface of a sphere with its unit 
vectors x,y,z pointing along the local 0, (j), r directions, 
respectively, as shown in Fig. [1] 

We solve the following set of MHD equations. The con- 
tinuity equation is given by 



Dlnp 
Dt 



(2) 



where D/Di = d/dt + U ■ "V denotes the Lagrangian 
derivative with respect to the local velocity of the gas U. 
Assuming an ideal gas, we express the pressure in terms of 
density, specific entropy s, and sound speed Cg, which, in 
turn, is a function of p and s. Thus the momentum equa- 
tion in a frame of reference rotating with angular velocity 
CI reads 



BU 



—— = - c, V h In p - 2r2 X (7 + .g^i 



J X B 

P 



1 



VV-J7 + 2S- Vlnp 



where J is the current density, B is the magnetic field, v is 
the constant kinematic viscosity, and S is the traceless rate- 
of-strain tensor. The sound speed is related to temperature 



by cj 



Cy)jT with Cp and c„ the specific heat at 



constant pressure and constant volume, respectively, and 
7 = Cp/cy is here fixed to 5/3. The induction equation is 
solved in terms of the magnetic vector potential A, such 
that ^ X A = B, hence 



dA 

'dt 



U X B + T]W^A, 



(4) 



where 77 denotes constant molecular magnetic diffusivity. 




Fig. 1. The Cartesian simulation domain with respect to 
spherical coordinates. 



Finally, we have for the entropy equation with temper- 
ature T and constant radiative (thermal) conductivity K 



pT 



Dt 



V • (KVT) + -qn^J^ + 2pvS'^ 



(5) 



where the temperature is related to the specific entropy by 



So + Cy In 



T/To 



{p/poY 



(6) 



We use the fully compressible Pencil CodeQ for all our 
calculations. 

For all quantities, periodic boundary conditions in the 
X and y directions are adopted. In the z direction we use 
the no-slip boundary condition for the velocity, the vertical 
field condition for the magnetic field, as a proxy for vacuum 
boundaries. We keep the temperature at the top and the 
(radiative) heat fiux at the bottom fixed. Their values were 
chosen to conform with the initial temperature profile of the 
(not magnetically modified) polytrope described below. 

2.1. Initial state 

The base state is a polytrope that is, p = Cp^ , with index 
m = l/(r— 1) = 3. The initial z profiles of density, pressure, 
temperature and entropy are given by. 



p, = po^Hz), p, = PQ<S>\z), T, = ro$(z), 
= So - Cy ln$(z). 



(7) 



where $ is a non-dimensional gravitational potential given 

by 

1 Qz 

1 + 



Hz) 



4 To(cp - Cy) 



|http://www.pencil-code.googlecode.com 
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with the reference point zq chosen to be at the bottom of 
the domain and the values at this point given by po, Po — 
c^qPo/7i Tq = C30/ (cp — Ci,)7 and sq. Here Cso is the reference 
sound speed to which we also refer to when calculating 
Mach numbers. 

As the adiabatic index is here mad = 1/(7^ 1) = 3/2, 
the subadiabaticity in the domain is very large, namely 
dlnT/dlnP- {dlnT/dln P)^d = -0.15. Thus, the initial 
stratification is highly stable to convection in the absence 
of any magnetic field, guaranteeing that turbulence is gen- 
erated solely by the buoyancy instability. 

The initial magnetic field is a horizontal layer of thick- 
ness Hb, where By has the profile 

Byo = So^B^ tanh (^^^) , (8) 

and the reference Alfven speed is defined by wao — 
Bq/ ^PafJ-o with /io being the vacuum permeability. If not 
indicated otherwise, the initial magnetic field strength is 
fixed to UAo/csO — 0.5. In order to satisfy the condi- 
tion ([T]) initially, we have to ensure Hg < Hp^zs), where 
Hp{z) — |Vlnp(z)|~^ is the local density scale height. 
When choosing zb — zi^ = O.SL^ this is satisfied for Hb < 
0.1Lz + 'iTo{cp~Cy)/3\gz\ which is surely true for the choice 
Hb = 0.05^2. 

Upon addition of a magnetic field, we have to modify the 
base state such that the density profile remains unchanged. 
In order to obey magnetostatic equilibrium, pressure and 
temperature are adjusted in the following way: 

P.^P.~^, T.^T.-I^ ' (9) 

ZfiQ ZfiQ Pi(Cp - Cy) 

The entropy is then re-calculated from Eq. ©. The initial 
velocity components Ux and Uy are specified such that it 
contains about 20 localized eddies in the plane z = zb 
with Mach numbers of about 10^^. Also the initial vertical 
velocity, Uz is Gaussian random noise with the same Mach 
number. The rms of the initial kinetic helicity, scaled with 
the product of initial rms velocity and vorticity, is denoted 

£K0, that is, £K0 = (V((W^-t^)'>Mm«W„„s)(0) - 4 X 

10-6. 

2.2. Control parameters, nondimensional quantities, and 
computational grid 

The problem posed by ([2]) through ([5]) is governed by five in- 
dependent dimensionless parameters, (i) the Prandtl num- 
ber Pr = vjxQi with the temperature conductivity xo — 
poCp^ (ii) the magnetic Prandtl number PrM = i^/?7, (iii) 
the "magnetic Taylor number" TaM = 2n'^Ly/rf, (iv) the 
rotational inclination (colatitude), 6, and (v) the normal- 
ized gravitational acceleration gzL^rf. In addition there 
are two independent parameters of the initial equilibrium 
(vi) the normalized pressure scale height at the bottom, 
Hp/Lz = c^g/^gzLz and (vii) the initial Lundquist num- 
ber, Luq = VAoHB/rj, based upon the thickness of the 
magnetic layer. In addition to this, we also have the non- 
dimensional sound speed, CsoLy/rj- this paper we shall 
keep the normalized pressure scale height and the sound 
speed fixed, while varying both Prandtl numbers, TaM, 
and Luq. The definitions as well as the values or ranges 
of the control parameters are summarized in Table [T] We 



Table 1. Non-dimensional control parameters character- 
izing the buoyancy instability. Note the definition of the 
modified plasma-beta /3 as the ratio of the total pressure 
Ptot = P + Pm to the magnetic pressure pu = ByQ/2po, 
because this quantity adopts a simple 1 /B'^ dependence on 
the magnetic field, cf Eq. ([9]) . Values of /3 refer to t = and 
the midplane of the magnetic sheet. 



Parameter 


Symbol 


Definition 


Value/Range 


norm, scale height 




Hp/Lz 


0.3 


norm, sound speed 




CsoLy/ri 


6 X 10'' 


Prandtl number 


Pr 




0.125 . . .4.0 


magnetic Prandtl no. 


PrM 


vjr] 


0.125 . . .4.0 


Roberts number 


Rb 


xoh 


0.25 ...1.0 


magnetic Taylor no. 


Taivi 




0...3.2 X 10^° 


rotational inclination 


e 




... 180 


(initial) Lundquist no. 


Luo 


vaoHb/v 


500 .. . 600 


(initial) modified 


Po 


(Ptot/PM)(ZB,0) 


1.04 . . . 3.22 


plasma-beta 







have also included in the same table two dependent pa- 
rameters namely the modified initial plasma-beta in the 
midplane of the magnetic layer and the Roberts number 
Rb = PrM/Pr-Xo/?y. 

The computational domain is defined by |a;| < 
\y\ < Ly/2, -Lz/4: < z< SLz/A, = Lz = Ly/3, thus its 
aspect ratio is 1:3:1. The results will be presented in non- 
dimensional form, velocity in units of the reference Alfven 
speed, uao: time in units of the corresponding Alfven travel 
time in the y direction, ^ao — Ly/v^o, and magnetic field 
in units of Bq or the rms value (J^ B^gdz/L^)^/^. 

It is instructive to look upon the relevant definitions 
of the fluid Reynolds number. Re, and the magnetic 
Reynolds number, Rcm, for this problem where the tur- 
bulence is driven solely by the instability of the magnetic 
layer. From first principles, the Re characterizes the ra- 
tio of the advective term {{U • 'VU)^)^^^ and the viscous 
term ((:/V^C/)^)^/^ in the Navier-Stokes equation, while 
Rejvi characterizes the ratio of ((V x {U x B)) )^/^ and 
{{rjV^B)"^)^/^ in the induction equation with the angular 
brackets representing volume averaging. Let us denote these 
ab initio definitions as "term-based" and refer to them by 
Re* and Re^. Note, that with the term-based definitions 
Rcm/Rc may well deviate from PrM- Alternatively, we can 
define a length scale Ljj — f/ims/STrWrms from the rms val- 
ues of velocity and vorticity and define the more conven- 
tional "length-based" Reynolds numbers Re — U^^sLu/v 
and Rcm = U^rasLu/r^. 

The calculations were carried out on equidistant grids 
with resolutions of either 64'^ or 128'^. For numerical testing 
we have also performed a few runs with 256'^ or 128^ x 256 
resolutions. 

2.3. The test- field method 

We now define mean magnetic and velocity fields, B and U, 
where overbars denote horizontal averaging. Fluctuations 
are defined correspondingly as h = B ~ B and u = U — 
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Table 2. List of runs of set B. The computational box is placed at colatitude 9 = 30°. Magnetic Taylor number 
TaM = 3.24 X 10^°, initial Lundquist number Luq = 500, initial plasma-beta /3q = 2.27 and resolution 128'^ throughout. 
Ma - Mach number, based on [/ims, - growth rate. Saturation reached at t^^'. For the mean EMF in the saturated 
stage global extrema of the dominating with respect to z and t are given. 



Run 


Pr 


PrM 


Ma 


UJl f AO 




10* X ST^'/vAoBo 


length-based 




term-based 














min 


max 


Re(Lc/) 


Re* 


ReM 


ReM/Re* 


B128a 


4.0 


4.0 


0.017 


15.6 


1.99 


-1.01 


2.34 


0.5 


0.4 


2.3 


5.8 


B128b 


1.0 


4.0 


0.036 


21.6 


1.42 


-3.39 


7.32 


0.9 


0.6 


2.8 


4.5 


B128c 


1.0 


1.0 


0.020 


13.2 


1.64 


-1.49 


3.03 


1.8 


1.4 


1.9 


1.4 


B128d 


0.25 


1.0 


0.038 


25.2 


1.27 


-4.02 


7.52 


2.9 


2.1 


2.8 


1.3 


B128e 


0.125 


0.5 


0.036 


24.0 


1.22 


-5.47 


6.19 


3.6 


3.3 


2.9 


0.9 


B128f 


0.125 


0.125 


0.043 


19.9 


1.54 


-3.50 


4.84 


8.2 


16.1 


3.1 


0.2 


B128g 


0.5 


0.5 


0.018 


19.2 


1.72 


-2.06 


3.69 


2.9 


2.5 


1.9 


0.8 


B128h 


0.5 


1.0 


0.032 


21.6 


1.67 


-3.94 


3.97 


1.7 


1.9 


3.2 


1.7 



U . Following the above convention, the induction equation 
may be horizontally averaged as. 



dB 

In 



= V X {U X B) +V X e + r]V^B, (10) 



where rj is the molecular magnetic diffusivity of the fluid 
(here assumed uniform), while S = u x b is the mean 
electromotive force. The essence of mean-field magneto- 
hydrodynamics is to provide an expression for 5 as a func- 
tion of the large scale magnetic field and its derivatives. 
Mathematically, 



£ = aB rjVB, 



(11) 



where a. and r] are called transport coefficients. Note that 
a much more general representation of £ is given by the 
convolution integral 



£{x,t) 



to 



G{x,x',t,t')B{x',t')d^x'dt' (12) 



with an appropriate tensorial kernel G. The aim of the test- 
field method is to provide an expression for G as a function 
of fluid properties. By subtracting the horizontally averaged 
equation from the real one, we obtain the following equation 
for the fluctuating magnetic field b. 



V X (U xbP" + UX BP' 



rjV^bP", (13) 



with, e^"? = ux hP'^ — The superscripts pq indicate that 
this equation is solved for suitably chosen test fields B'P'^ 
with p, q = 1 , 2 if a and "q are assumed to be 2 x 2 matrices. 
This is the equation invoked by the test-field method for 
calculating the tensors a. and -q. The test- field suite of the 
Pencil Code has the provision for using either harmonic 
test fields i.e., 



= (cosfcz,0,0), B^^ = (O,cosA:0,O), 

21 22 

B = (sinA;z,0,0), B = (0, sin fcz, 0), 



(14) 



or linear test fields i.e. 



11 12 

B =(1,0,0), B =(0,1,0), 
:b'' = (z,0,0), #' = (0,z,0). 



(15) 



When it comes to applying the test-field method, an as- 
pect not discussed up to now is the intrinsic inhomogeneity 
of the flow both due to stratification and the background 
magnetic field itself. Within kinematics, that is without the 
background field, no specific complication is connected to 
this as a and r] emerge straightforwardly from the sta- 
tionary version of Equation (1121) in a shape expressing in- 
homogeneity, that is, a.{x,x'), ri{x,x') or, equivalently, 
a{x,x — x'), r]{x,x — x'). Performing a Fourier transform 
with respect to their second argument, we arrive at a{x, k) 
and f){x,k). In our case, harmonic test fields with differ- 
ent wavenumbers k in the z direction can be employed to 
obtain a{z, k) and r)(z, k). 

In the nonlinear situation, the Green's function ap- 
proach remains valid if £ is considered as a functional of 
U and B which is then linear and homogeneous in the lat- 
ter. However, we have to label G by the B actually acting 
upon U, that is, G{x, x'; B), and can thus only make state- 
ments about the transport tensors for just the particular 
B at hand. Hence, the tensors have to be labelled likewise: 
a{z, fc; B), 'i){z, k; B). As our initial mean magnetic field is 
in the y direction, the instability will generate a Bx and we 
are mainly interested in the coefficients a2i, 022, ?72i and 
7722 with rank-2 tensor components rjij = —rjikse-jks- 



3. Results 

3.1. Nature of the instability 

To start with we have performed a number of runs with 
different values of Pr and PrM, but all other dimensionless 
parameters held fixed, see Tabic [21 In particular, we have 
used a value of TaM = 3.24 x lO^'^ for the magnetic Taylor 
number and Luq = 500 for the initial Lundquist number. 
Table [2] shows the Reynolds numbers according to the two 
alternative definitions provided in Section 12.21 Note that 
with the exception of the run B128f, Re from the "length- 
based" and the "term-based" definitions are in agreement. 
Also the ratio RcM/Re* from the term-based definitions 
approaches PrM reasonably. 

We first show the temporal evolution of the magnetic 
field for a few representative cases in Fig. [51 In all of 
them, we can clearly distinguish a first stage of exponential 
growth, from a subsequent saturation phase. The x and z 
components of the magnetic field are generated at the ex- 
pense of its y component. Although there exists a persistent 
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energy source in the form of a constant heat flux into the 
domain, the final saturated stage always undergoes a slow 
decay. This decay is most clearly visible in By. Thus the 
instability is not able to maintain a dynamo on its own. 

We suppose that the magnetic layer formed by B^, 
though having a vertical scale suited to maintain the in- 
stability, is eventually not strong enough to take over the 
role of the initial magnetic layer. Let us first discuss the 
initial linear stage of the instability. 

3.2. Linear stage. 

At first we verify that the instability is indeed driven by 
magnetic buoyancy. As the coefficients in Eqs. ([l])-® are 
constant, the initial state ([7]) depends only on z, and the 
boundary conditions in the x and y directions are periodic, 
all eigensolutions ij) E {p, u, b, s} of the linearized problem 
must have the form 

^P{x, t) = iP{z) e' [2vr(mx/L, + ny/Ly) - u;t] ^ (^g) 

where m and n are integers and uj — cjr + iioi. 
Corresponding dispersion relations uj{m,n) have been es- 
tablished by applying perturbations of the forn i (1161) with 
a variational principle in the non-rotating case (jFai ri2nnih 
and with the set of linearized anelastic MHD equations in 
magnetostrophi c approximation at finite angular velocity 
(ISchmitti . 1198511 . The former case allows both oscillatory 
and non-oscillatory unstable modes, although in lFanI (|2001[ ) 
only non-oscillatory modes are reported. In the latter case, 
however, all unstable modes turn out to be oscillatory with 
the ratio ujr/u ji decreasing wi th latitude. Note that the ana- 
lytic results of lSchmit3 ()1985[ ) are limited in their predictive 
power by the fact that the variables are not subjected to 
our specific boundary conditions and that the analysis is 
performed locally. 

For the runs in Table [5] we find that in the early expo- 
nential growth phase m = 8 and n — 1 throughout as seen 
in Fig. 13] which shows a typical velocity pattern at a time 
during the linear stage. 

This is consistent with the findings of lFanI (j200lt ) where 
the fastest growing mode had always the smallest possible 
(non-vanishing) wavenumber in the direction of the field 
whereas the wavenumber perpendicular to the field was 
high. According to the terminology of Fan we may qualify 
our eigenmodes as undular as they change periodically in 
the direction of the magnetic background field. In our case, 
there seems to be some mixing with lower m modes since 
the growing perturbations do not appear to be perfectly si- 
nusoidal. While the growth rates presented in Table[2]could 
be easily identified from the averaged quantities shown in 
Fig. [21 it was difficult to access the oscillation frequencies. 
This is because they are small compared to the growth rates 
and saturation sets in too early to allow for the observation 
of a complete oscillation period. Nevertheless, some indica- 
tions for temporal variations in the eigenmode geometries 
have been found. 

Generally, we observe an increase of the growth rate 
with increasing magnetic Prandtl number, but a decrease 
with increasing Prandtl number. We find that the growth 
rate increases with the Roberts number as shown in Fig. 21 
This means that increasing efficiency of heat conduction 
in comparison to magnetic diffusion destabilizes the sub- 
adiabatic stratification in the system in agreement with the 
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Fig. 2. Time evolution of the runs in Table [H 
Upper panel: rms values of velocity and generated 
magnetic field components B^, B^ of run B128c 
(Pr = PrM = 1) scaled by (u^g) and (B^o) respectively. 
Note the clear exponential growth until t w 1.4iAo- Fast 
oscillations in ([/^) until t ~ ^ao indicate g-modes origi- 
nating from the initial velocity perturbation. Middle panel: 
rms values of generated magnetic field components for 
different runs. For legend see lower panel. Prandtl numbers 
indicated as (Pr,PrM). Lower panel: rms value of By. 



destab ilizing effect of thermal diffusion studied bv lAchesonI 
(IT979h . 
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3.3. Dependence on initial magnetic field and rotation 

Another piece of evidence for the magnetic character of 
the instability is its dependence on the initial magnetic 
field strength. From Fig. [5] we see a clear increase of the 
growth rate and saturation level with decreasing /Sq, that 
is, increasing Luq, while keeping the rotation rate fixed at 
TaM = 3.24 X 10^". Schmitt (2000) predicted a growth rate 
oc v^q/H for finite rotation, in the magnetostrophic approx- 
imation, inversely proportional to C = I^oVTom- 

Next we keep /3q constant at 2.27 and decrease ( grad- 
ually from 4.09 x 10^ to 0. Inspecting Fig. [SJ we find that 
growth rate and saturation level of (i?^ + B'^) increase 
monotonically and reach their maxima at ^ = (fi = 0) 
while the saturation time is decreasing. The impeding ef- 
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Fig. 3. Top: Velocity components Uy (in color), Ux and 
Uz (vectors) in the plane y = 0. Bottom: Ux (in color), Uy 
and Uz (vectors) in the plane x = 0. Both during the linear 
evolution phase of the run B128a. 
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Fig. 5. Dependence of the instability on initial magnetic 
field strength, expressed by /3q and rotation, expressed by 

— 1 /2 

= P^TaM ■ Upper panel: rms value of the generated 
magnetic field components (i?^ + Bf). Lower panel: {By). 

Legend shows (/3q,10^^C)- Pi' = Pi'm = 1, colatitude 9 = 
30°, resolution 64^ throughout. Note that the normalization 
time ^Ao is not the same for all curves, but proportional to 
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Fig. 4. Dependence of growth rate wi on the inverse 
Roberts number derived from the runs in Table [2] Solid 
line: best linear fit. Size of circles codes for the value of Re 
(length-based, see Tab. [5]). 



feet of rotation onto the instability at large ft is plausible 
in view of the Taylor-Proudman theorem because the un- 
stable eigenmodes do show pronounced z gradients in U, 
see Fig. il 



3.4. Saturated stage. 

At later time the instability reaches saturation, character- 
ized by turbulent magnetic, velocity, density and tempera- 
ture fields, that decay slowly thereafter. However, in most 
of the analysis below, this decay will be ignored and the 
turbulence approximately statistically stationary. The tur- 
bulence is necessarily both inhomogeneous and anisotropic 
and we shall further show that it is also helical. Under such 
conditions we expect the emergence of a mean electromo- 
tive force. Indeed magnetic fields perpendicular to the ini- 
tial magnetic layer are produced having non-vanishing hor- 
izontal averages. 
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Pr = 0.125, PrM = 0.125, Re = 8.2 




Pr = 0.125, PiM = 0.5, Re = 3.6 





Fig. 6. Left: volume rendering of the By = Q.IBq isosurface. Right: field lines, colored according to the value of By for 
runs B128f (top) and B128e (bottom) a.t t = 2tAo (saturated stage). 



In order to give a better idea of the 3D geometry of 
the magnetic field we provide in Fig. [5] a volume rendering 
of By at a time after t'''^*- for the runs B128e and B128f 
(see Table (H) which differ only in their magnetic Prandtl 
numbers. Notice how the magnetic layer breaks into fiux 
tubes - similar to what is see n in Fig. 3 of Fan (2001) and 
also in [Matthews et all (|1995[ ). The difference between the 
two cases is most striking in the nature of corrugation in the 
surface shown. We attribute the difference to larger twist in 
the rising tubular structures in the run B128e compared to 
B128f which becomes clearly visible in the field line pictures 
also depicted in Fig. [5] (right). 



Figure [7] demonstrates the breakup of the magnetic layer 
into tubular structures of concentrated magnetic field which 
are also regions of low density, hence rising. Notice also 
the high density regions just above and below these tubu- 



lar structures. They show a significantly lower temperature 
than their surroundings (Fig. [71 bottom). 



Considering the solar convection zone it is suggestive to 
ask to what extent the flux tubes are twisted, as their abil- 
ity to rise over a large distance depends crucially on this 
property. For a quantitative measurement we utilize the di- 
mensionless parameter ej — {J ■ S)/ Jinis^rms, the relative 
current helicity, essentially measuring the overall degree of 
alignment between B and J. Here, angular brackets de- 
note volume averages. A corresponding localized quantity 
is ej{z) = J ■ B/ JimsBrms. Figure [5] shows sj (filled con- 
tours) as well as By in the plane y ^ for run B128e. Notice 
that the contours are bend leftward because of the Coriolis 
force with Q ■ z > 0. The contour plots of By in this figure 
also show the formation of rising tubular structures from 
the magnetic layer. 
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Fig. 7. Relative density perturbation, 8p/ pi (top) and 
relative temperature perturbation 5T/Ti (bottom), with 
5p = p — Pi, ST = T — Ti and Pi{z), Ti(z) taken from 
Eq. d?]), in the plane y = at t = t""^*- ^ 2tAo for the run 
B128h. Both plots overlaid with contours of By (solid lines). 

In Fig. |9]we show the dependence of ej on Rb and the 
profiles ej{z) for some selected runs. Although the total 
helicity reaches only values of a few percent, its localized 
counterpart is as strong as 30% near to the initial location 
of the magnetic sheet. The clear dependence of ej on Rb is 
in contrast to the only weak dependences on PrM and Pr 
individually. This is an important result from this section. 
Our conjecture is that, at large Rb, this magnetic buoyancy 
instability may play an important role in the formation of 
twisted flux tubes in the Sun, where Rb ^ 1 is expected. 

To demonstrate the emergence of a mean magnetic field 
we present in Fig. [TUltime-depth plots of By and B^ for the 
run B128g (note that B^ — 0). There, t « 1.6 iAo marks the 
end of the exponential growth phase after which a strong 
growth of Bx, obviously at the expense of By, sets in. B^ 
reaches its maximum around t « 3tAo and is then subject 
to the overall decay. Note the strong vertical concentration 
of Bx, approximately antisymmetric about the midplane of 
the magnetic sheet. 



3.5. Calculation of turbulent transport coefficients 

The turbulence resulting from the buoyancy instability gen- 
erates a mean magnetic field component B^ from an initial 
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Fig. 8. Scaled current helicity J ■ B/Jrms-Brms for run 
B128d (top) at i = t^'^^ ~ 2tAo- Arrows show and v^. 
By/Byo (bottom) for the same run. Arrows show B^ and 
Bz. Both panels show the plane y = 0. 
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Fig. 9. Top: Dependence of the total relative current helic- 
ity sj on inverse Roberts number. Size of circles codes for 
value of Rcm- Bottom: Dependence of J ■ J5/ Ji-nis^rms on z 
after saturation for runs B128e (solid, Rb"^ = 0.25), B128h 
(dashed, Rb^^ 0.5), and B128c (dotted, Rb"^ = 1). 

By which is also modified compared to its initial shape 



(see Fig. [TU]). It is then natural to employ the technique 
known as the quasi-kinematic test-field method to calcu- 
late transport coefficients like the a. and tensors which 
describe this process. So far, test-fields have mostly been 
used in situations where a hydrodynamic background was 
already present in absence of the mean magnetic field (see. 
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e.g. iBrandenburg et all l2008al lblR). Here, in contrast, the 
(magnetohydrodynamic) turbulence results entirely from 
the instability of a pre-existing mean magnetic field, Byo{z). 
In other words, our simulations do not posses a kinematic 
stage in which the influence of B would be negligible. One 
might worry that in s uch a situation the quas i-kinematic 
test-field method fails (jCourvoisier et al.L 1201011 . However, 
Eq. (|13l) continues to be valid and hence all conclusions 
drawn from it, because the decisive applicability criterion 
is whether or not there exists hydromagnetic turbulence in 
the absence of the mean magnetic field. This is not the case 
here, so the method should be applicable. The only pecu- 
liarity occurring is the fact that all components of ol and 
77 vanish for < B^ms < -Bthreshoid, because fluctuating 
velocity and magnetic flelds develop only after the insta- 
bility has set in. Another aspect not considered in most 
previous test-fleld studies is the strong intrinsic inhomo- 
geneity of the turbulence not only as a consequence of the 
strong z dependence of B, but also due to the stratifled 
density background. Thus the transport coefficients need 
to be determined as z dependent quantities. We shall next 
demonstrate that the test-field method still works reason- 
ably well in this regime. Note that to calculate the transport 
coefficients in addition to the usual MHD equations four 
additional evolution equations of the form (1131) for four in- 
dependent test-fields have to be solved. Hence the test-field 
runs are computationally almost thrice as expensive. We 
have thus reduced resolution to 64^ grid points for all these 
runs. 



0.15 
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(b) BxlO^ 




Fig. 10. Time-depth diagram for Bx,y normalized on Bq 
for run B128g in Table[2](2; extent of the box clipped). Note 
the difference of two orders of magnitude between By and 
B^. 



3.5.1. Reconstruction of the mean EMF 

To validate the test-field method we first confirm that the 
quantity £, taken directly from the DNS, can be reproduced 
by employing the relation (jll[) between £ and B with the 
tensors at and J7 determined using the quasi-kinematic test- 
field method. In mathematical terms. 



£f{z-B) 



with 



E 

k 



[tl^{z, k- B)B]{k) + iClj{z, k- B)B][k)\ (17) 



ICIj = aij{z, k"; B) cos{k''z) - f]ii{z, k"] B^ij^k" sin(fc''z), 
^] = ajj(2;,fc^B)sin(Fz) + Tiu{z,k'';B)eij3k'' cos{k''z), 
2 



Bj{k) = — / Bj(z)cos(fc^£)dz, 

Lz Jz 

B]{k) = ^ / Bj(z)sin(Fz)dz, 

Lz Jz 



(18) 



Lz Lz 



Z- Zq- 



2 ' 



where the superscript R indicates reconstruction. Here, the 
boundary condition for B gives rise to the selection of dis- 
crete cosine and sine modes with wavenumbers k'^ and k^ ^ 
respectively. The additional argument B is to indicate that 
the kernels /C°'**, as well as the tensors a. and 77, are valid 
just for that mean field B which is present in the main run. 
As a consequence, the reconstruction of the mean EMF can 
be successful only when employing exactly this B in (jl8p . 
That is, the mean field representation of the turbulence by 
a. and r\ has, at this level, merely descriptive rather than 
predictive potential. 

Let us denote 5^ as the reconstructed EMF according 
to Eq. (IT7I) truncated at k' < fcmax, with fc^^^ — 2k'Ti/Lz. 
Here k' can now take both integer and half-integer values 
where the integer (half-integer) values of k' correspond to 
the family of sine (cosine) modes in Eq. (ITSt . An initial 
estimate of fcmax required for a reasonable reconstruction 
of £ was obtained frorn_the power spectra of both B^ and 
By. It turned out that B^ has significant spectral power up 
until fc' — 16, whereas for By the power spectra has levelled 
off already at fc' = 8. The components of the tensors ol and 
77 also show rather different spectral behavior, both in the 
midplane of the magnetic layer and near the midplane of 
the box as seen in Fig. 1111 From the figures it is evident 
that in most cases the spectra can be reasonably truncated 
at fc' = 8 with the exceptions of ai2 and a2i- Note that 
the values for fc' = are not relevant here as, due to the 
boundary conditions, B does not posses_s a fc' = contri- 
bution. The result of the assembly of £^ from (fT7|) with 
(fT5)) . is presented in Fig. [T^ middle column. From simple 
visual inspection we find it to be a faithful reproduction of 
£ from the DNS shown in the left column. Clearly, a naive 
application of the test-field procedure with harmonic test 
fields with only the lowest fc' — 0.5 results in an inadequate 
description as shown in the right column. We define two 
measures for the quality of the mean EMF reconstruction 
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Fig. 11. Dependence of a and rj on the test-field wavenumber k' in the midplane of the magnetic layer (z = zb) and 
near the midplane of the box {z = 0.21Lz) for run TF30+ of Table. [31 Integer and half-integer values of k' belong to sine 
and cosine modes in B, respectively. Note that fc' = refers to constant and linear test fields and that the coefficients 
for that value do not enter the S - B relation for the given setup. 



namely xt' ^^'^ correlation coefficient rk' defined as 

Xl = r„ = , (19) 

^ {Sly}. AS^l}.,t 

where the subscript "z,t" denotes that the averaging has 
been carried out over the vertical coordinate z as well as 
over the temporal range 1.2<ao < ^ < 3.4tAo- The rela- 
tive error of the reconstruction, x|, , and the correlation 
coefficient, rk', are plotted in Fig. [T3] as a function of the 
truncation wavenumber fcmax- The xt' reach a minimum 
value and level off around k' — 8 for both and £y. This 
implies that including higher harmonic test fields beyond 
k' = 7 does not improve the reconstructed EMF. We spec- 
ulate that the reason behind t his discrepancy is that we 
have neglected memory effects (jHubbard fc Brandenburel 
l2009tl in the turbulent transport coefficients. This can be 
particularly important in the present situation as we are 
obviously not in a statistically stationary regime. Similarly 
rk' for £^ {£^) converges to a value of 0.98 (0.93) at fc' = 4 



(8). It is important to note that even though the tensor 
components ai2 and a2i do not converge with increasing 
k', the reconstructed EMFs do. Also calculating transport 
coefficients for k' > 8 does not improve the reconstruction 
any further. This is probably because we do not sufficiently 
resolve wavenumber scales larger than 10 in the domain 
with a grid resolution of only 64"^ . 

3.5.2. Dependence of the transport tensors on inclination 

From the point of view of the solar dynamo it is important 
to look at a and rj as functions of the rotational inclina- 
tion 6 or latitude A, with a focus on symmetry properties 
with respect to A = 0, which is the solar equator. Moving 
from the northern hemisphere at A to the southern at —A, 
that is changing 6 to tt ~ 9, but keeping all other problem 
parameters constant, is equivalent to inverting the sign of 
^Iz . As the same can be accomplished by reflecting the cor- 
responding rigid rotation about the plane x = 0, we might 
construct the solution (p, U, B, s) of ([1])-® for —A simply 
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Fig. 12. Reconstruction of the mean EMF for the run 
TF30+ using a and 77 from the test-field method. Top: 
£x{z,t), bottom: £y{z,t), both scaled by IO^^vaoBq- Left: 
directly from u x b. Middle: Reconstruction using all con- 
tributions k' = 0.5, 1, 1.5, ...16 in (fT7|) . Right: Same as be- 
fore, but using only the k' = 0.5 contribution. 
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Fig. 13. Quality of the EMF reconstruction as a function of 
the truncation wavenumber fcmax^ Xfe' (left) and correlation 
rk' (right) calculated for £^ (solid) and £^ (dashed) using 
Eq. (Ull). 



by reflecting it properly about the same plane. Under this 
reflection polar vectors like velocity transform as, 

{f/„ Uy, C/J {x, y, z) ^ {-C/„ Uy, C/J {~x, y, z), (20) 

and axial vectors like the magnetic fleld as 

{B,,By,B,}{x,y,z) ^ {B,,^By,B,}{-x,y,z) (21) 

(Note that the gravitational acceleration is invariant un- 
der this reflection.) Hence, for the initial magnetic field, 
Byo{z), the transition to —A requires only a sign inver- 
sion. But, since the induction equation is linear in B, and 



Lorcntz force as well as Ohmic dissipation are quadratic, 
inverting the sign of ByQ{z) would just transform the so- 
lution {p,U, B, s} to {p,U,—B,s}, that is, would leave 
the turbulence essentially unchanged and can be omitted. 
Moreover, as the transport coefficients, expressing correla- 
tion properties of the turbulent velocity it, are functions of z 
only the reflection operation can hardly change their magni- 
tudes. With respect to possible sign inversions we consider, 
that £ and J, being polar vectors, invert the sign of their x 
components under reflection, but keep their y components 
unchanged. The axial vector B behaves just the opposite 
way. Thus, we have an — >■ —an for z = 1, 2 (no summation) 
and aij — )■ a^- for i ^ j, whereas rju — > -qu for i = 1,2 
and rjij — )■ —i^ij for i ^ j when moving from A to —A. 
Consequently, it appears that the results for the southern 
hemisphere can be derived from those for the northern by 
simple operations. Strictly speaking however, this is only 
true when the initial condition for U is also reflected upon 
the transition from A to —A. From a naive point of view we 
might suppose that omitting this reflection can hardly be of 
any importance, because we use random initial condition. 
But this we have found not to be true. We note further that 
once the initial condition is reflected too the symmetry is 
restored. 



According to the resuhs of ISchmittI ()2000t) we expect 
a decrease in the intensity of the instability with increas- 
ing inclination of the rotation axis. This can be explained 
by the buoyant nature of the turbulence, for which verti- 
cal motions are essential. At the poles, the effect of the 
Coriolis force on vertical motions is weakest, whereas they 
are strongly deflected at the equator. Figure [Til indeed con- 
firms, that the growth rates decrease continuously when 
changing 6 from 0° towards 90°. In Fig. [TS] we show the 
variation of the mean magnetic field and the corresponding 
mean EMF with latitude and z at a time during the satu- 
rated stage. A pecularity in this figure is that £y and con- 
sequently Bx are non-zero at the equator where we would 
expect these quantities to vanish. This is an example of 
spontaneous symmetry breaking and can be explained by a 
mean field dynamo operating at the equator. This dynamo 
generates B^ whose sign is determined by the random ini- 
tial conditions. A detailed discussion of this issue will be 
provided in a forthcoming paper. In the rest of the paper 
we anti-symmetrize Bx and £y, while symmetrize By and 
£x about the equator (see Fig. [121). This is done by including 
the results from runs with two different initial conditions 
for velocity, one being the mirror reflection of the other ac- 
cording to Eq. d^ni- In particular at the equator {9 — 90°), 
the two initial conditions give rise to a Bx with exactly 
the same magnitude but differing in sign. Thus averaging 
the Bx from the two runs gives a zero Bx at the equator. 
We perform the same operation for the turbulent transport 
coefficients calculated from the QKTF method. The trans- 
port coefficients calculated from only the k' — 0.5 test fields 
belonging to the family of cosine modes are presented in 
Fig-dZl It can be immediately seen from these plots that the 
instability becomes more effective with increasing (north- 
ern or southern) latitude. Corresponding runs performed 
with linear test- fields are compiled in Table [31 We observe 
that the turbulent transport coefficients increase in modu- 
lus when moving towards the poles, but are, with the only 
exception of 0121, significantly reduced close to the equator. 
Obviously, the transport coefficients respond directly to the 
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Table 3. List of runs from set TF. Pr = 4.0, PrM = 4.0, 
TaM = 3.24 X 10^", except for last three runs with TaM — 0. 
^0 = 2.27 except in last two runs, TFOOl and TFOOm, where 
/3o is 3.22 and 1.03, respectively. Resolution 64"^ throughout. 
Saturation is reached at t^^^ . Global extrema of the dom- 
inating Ey with respect to z and t are given. Re ~ 4.8 
throughout. The arrows in the f2 column indicate the sign 

of riz. 
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Fig. 14. Dependence of the instability on rotational incli- 
nation Q in terms of rms value of generated field components 
(Bl + Bl) for the runs TFO+, TF30+, TF60-h and TF90-f 
in Table. El 



inhibition of the vertical motions by the Coriolis force when 
moving towards the equator. 

The a tensor can be decomposed in symmetric and an- 
tisymmetric parts. The latter represents a turbulent pump- 
ing velocity 7, and gives rise to the term 7 x B in the 
mean EMF. By virtue of the horizontal averaging of the 
magnetic field, = 0. Hence, the only relevant compo- 
nent of pumping is 7^ which is defined by (a2i ~ Q!i2)/2. 
Analytical results indicate that in a wide range of situa- 
tions, the turbulent pumping is directed away from the re- 
gion of strong turbulence ("turbulent diamagnetism" , see 
iKrause fc Radler, 1980). From Fig.[Tl] we see that the com- 
ponents a2i and ai2 do not converge to zero with increas- 
ing k' . In fact ai2 changes sign at fc' = 2 and 021 does 
so at fc' = 8. Consequently 7^ determined from harmonic 
test fields with fc' — 0.5 and fc' > 4 should have opposite 
signs as confirmed by Fig. [TS] Physically, this means that 
magnetic fields formed on the scale of Lz will be pumped 
away from the initial magnetic layer while those on the 



Fig. 15. Latitudinal dependence of Bx.y{z), £x,y{z) aver- 
aged between t = t^^*" and t^^^+t ao- All problem parameters 
except 9 held fixed at the values of run TFO-I- in Table [31 
Dashed line: initial position of the magnetic layer. 
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Fig. 16. Same as Fig. [151 but now symmetrized about the 
equator. 



scale of the magnetic layer, Hb shall be pumped into the 
layer, the latter being contrary to the standard concept 
of "turbulent diamagnetism". It is thus difficult to com- 
ment on the transport of the total B by 7. Only if the 
pumping were oriented away from the magnetic layer for 
all the wavenumbers of the dominating constituents in B it 
would lead to a broadening of the initial layer i.e., a reduc- 
tion of dBy/dz and would hence inhibit the instability. A 
similar dependence of turbulent pumping on wavenumber 
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Fig. 17. Dependence of a. and r;, averaged between t = t^^^ 
and t'*'^* -I- on latitude, A = 90° - 9, and z, calculated 
using test fields with k' — 0.5. All problem parameters ex- 
cept 9 held fixed at the values of run TFO-I- in Table 2. 
a scaled by i/rms, J? scaled by the molecular diffusivity 77. 
Dashed line: initial position of the magnetic layer. 
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Fig. 18. Pumping velocity 7^ — (a2i — ai2)/2 scaled with 
Urms, as a function of latitude, A for k' ^ 0.5 (left panel) 
and k' ~ 4.5 (right panel). 



has been found bv lKapvla. Korpi fc BrandenbursJ (|2009f ) in 
DNS of convection. With regard to to the saturation of the 



magnetic buoyancy instability, a strong turbulent magnetic 
diffusion given by 7722 (see Fig. I14p is likely to be more im- 
portant. At the poles this quantity is as large as 40 times 
the molecular value of rj. 



4. Conclusions 

We have studied in detail the generation of the a effect due 
to the buoyancy instability of a toroidal magnetic layer in 
a stratified atmosphere by using direct numerical simula- 
tions. We find that both the magnetic energy and the cur- 
rent helicity in the system increase monotonically with the 
ratio of thermal conductivity to magnetic diffusivity, the 
Roberts numb er Rb (F i g.Bt. This agrees wi t h ear lier ana- 
lytical work of iGilmanI (llQTOT) and lAchesonI (| 19 791 ) as well 
as numerical work of ' iSilvers et al.l (I2009D which find that 
efficient thermal diffusion or heat exchange can destabilize 
a stable stratification. The dependence of twist on Rb is an 
important result since the buoyancy instability would pro- 
duce twisted flux tubes from a m agnetic layer, if it e xisted 
in the overshoot layer of the Sun. 'Vasil fc B rummel (j2008|) 
also reported the formation of twisted fiux tubes from a 
horizontal magnetic layer produced, but in their case it is 
due to the action of shear on a weak vertical magnetic field. 
We further find that the growth rate of the buoyancy in- 
stability is reduced in presence of rotation compared to the 
case with i7 = 0. 

We have run our simulations only until the time taken 
by the initial magnetic layer to break up due to the buoy- 
ancy instability. In absence of any other forcing such as a 
strong shear, the buoyancy instability cannot usually sus- 
tain itself past the break-up phase since the vertical gradi- 
ent of the magnetic energy in the layer becomes compara- 
ble to the stratification due to magnetic diffusion. We may 
say that strong shear is not imperative to the production 
of tubular structures from the toroidal magnetic layer but 
will play a key role in keeping the layer from breaking up. 
It may also be possible that turbulent pumping arrests the 
decay of such a magnetic layer in the actual overshoot re- 
gion. However, it is not yet clear if such a layer exists and 
is subject to the buoyancy instability in the real Sun. 

We have 'measured' the turbulent transport coeffi- 
cients using the technique of the quasi-kinematic test- 
field method. In order to prove that the a and rj ten- 
sors obtained from this method are reasonably accurate, we 
show the agreement between £ = u' x b' and the ansatz 
£ = aB — r/J using harmonic test fields with wavenum- 
bers < fc' < 16. Here we have illustrated a technique 
of judging the reliability of transport coefficients obtained 
from the test- field method. We find that, even in presence 
of magnetically driven turbulence, a and rj obtained from 
the quasi-kinematic test-field method provide a reasonably 
accurate description of the turbulent EMF. This is an im- 
portant outcome of our study. 

We find that determined using a harmonic test field 
with the lowest wavenumber that fits in vertical extent of 
the box already comprises a considerable part of the total 
EMF. Hence we can use QKTF to calculate the turbulent 
coefficients at finite J7 as a function of latitude using har- 
monic test fields with this wavenumber. The component 
Q22 contributes to the generation of from the strong 



initial field By in the layer. The off-diagonal components 
contribute to a vertical turbulent pumping velocity directed 
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away from the region of turbulence surrounding the mag- 
netic layer. The influence of this component systematically 
expands along z with increasing latitude and somewhat 
agrees with the result in Brandenburg & Schmitt (1998). 
The agreement is not complete since the a.2i{z^ff) is inho- 
mogeneous with respect to z and can have sign changes 
along 6*, e.g., at zjL^ = 0.4 in Fig. [TTJi. We find that ah 
transport coefficients except a2i increase with latitude and 
are significantly reduced near the equator due to the sup- 
pressing effect of the Coriolis force on the instability. 

For the first time the turbulent magnetic diffusivity 
given by the diagonal components of r\ has been computed, 
as shown in Fig. 1171 In particular, near the magnetic layer, 
the diagonal component 7722 is 25 times larger than the 
molecular value 77. The buoyancy driven instability has the 
property that the a as measured by the growth rate of the 
instability increases with the magnitude of the magnetic 
field in the horizontal layer (compare solid and dashed lines 
Fig. [5]). This property makes it an attractive candidate for 
solar dynamo models, unlike the a generated due to he- 
lical turbulence which gets quenched for strong magnetic 
fields. The increase of a and r\ with B is a remarkable re- 
sult an d supports similar suggestions bv lBrandenburg et al.l 
(|1998D that, if turbulent transport coefficients are caused 
by fiows that are magnetically driven like here or, e.g., 
in Balbus-Hawley instabilities, then both a and 77 may in- 
crease with the magnetic field strength. This trend is some- 
times referred to as 'anti-quenching' and may be needed 
to support the observational relation between the ratio of 
dynamo cycle to rotation frequencies, Wcyc/i^ and Rossby 
number inverse, Ro~^ for stellar d ata (jBrandenburg et all . 
Il998t ISaar fc BrandenbureL I1999D . Note finally that mod- 
elling the a as a function of space and the mean magnetic 
field to use in a mean field dynamo model is a very difficult 
proposition that needs to be postponed to future work. 
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